A classification of scalar field potentials with cosmological scaling solutions 
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An attractive method of obtaining an effective cosmological constant at the present epoch is through 
the potential energy of a scalar field. Considering models with a perfect fluid and a scalar field, 
we classify all potentials for which the scalar field energy density scales as a power-law of the scale 
factor when the perfect fluid density dominates. There are three possibilities. The first two are well 
known; the much-investigated exponential potentials have the scalar field mimicking the evolution 
of the perfect fluid, while for negative power-laws, introduced by Ratra and Peebles, the scalar field 
density grows relative to that of the fluid. The third possibility is a new one, where the potential 
is a positive power-law and the scalar field energy density decays relative to the perfect fluid. We 
provide a complete analysis of exact solutions and their stability properties, and investigate a range 
of possible cosmological applications. 



PACS numbers: 98.80.Cq 



I. INTRODUCTION 
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The evidence in favor of a cosmological constant, or 
something very much like it, playing a significant dy- 
namical role in our present Universe is becoming over- 
whelming. Most prominent have been the recent mea- 
surements of the apparent magnitude— redshift relation- 
ship using Type la supernovae [fjj , but other factors such 
as the consistently low measurements of the matter den- 
sity, including the baryon fraction in galaxy clusters 
have also been pointing in that direction. While many of 
these latter measurements are insensitive to the presence 
of a cosmological constant, there is some observational 
motivation for a flat Universe from the favored location 
of the first acoustic peak of the microwave background 
anisotropics, and some theoretical motivation from a de- 
sire to utilize the simplest models of cosmological infla- 
tion as the source of density perturbations. In combina- 
tion, these favor a present cosmological constant (in units 
of the critical density) of VL\ ~ 0.7. 

Since a genuine cosmological constant requires extreme 
fine-tuning in order to have only begun to dominate re- 
cently, it is extremely tempting to model the cosmolog- 
ical constant as an effective one. As the supernovae ob- 
servations are requiring an accelerating Universe, which 
is precisely the definition of inflation, the minimal ap- 
proach is to assume that the same mechanism drives in- 
flation now as is presumed in the early Universe, namely 
the potential energy of a scalar field. Three possibilities 
present themselves. The field could be at an absolute 
minimum of non-zero potential energy It could be in a 
metastable false vacuum, tunneling at some later stage 
into the true vacuum and perhaps even reheating. Or it 
could be slowly rolling down a potential, as in the chaotic 
inflation models favored for the early Universe. 



The first two of these possibilities are dynamically in- 
distinguishable from a true cosmological constant, and so 
we shall concentrate on the third, which is often called 
'quintessence'. As stressed in a recent paper by Zlatev 
et al. H|, a rolling scalar field offers an opportunity to 
address another mystery, that of why the cosmological 
constant took so long to become dominant. If, for exam- 
ple, the scalar field behaves in such a way as to remain 
insignificant during the radiation domination era, per- 
haps it can be 'triggered' in some way to begin to grow 
in the matter era and come to dominate only in the re- 
cent past. Solutions where the scalar field energy density 
follows that of radiation or matter have been called "scal- 
ing solutions" , and more recently "trackers" , and several 
examples have been described in the literature. 

For the purposes of this paper, we will define a "scaling 
solution" as one in which the scalar field energy density, 
p$ scales exactly as a power of the scale factor, p$ oc R~ n 
when the dominant component has an energy density 
which scales as a (possibly different) power: p oc R~ m . 
An equivalent, and perhaps more fundamental, defini- 
tion is that the scalar field kinetic and potential energies 
maintain a fixed ratio. We will use the term "tracker 
solution" to refer to the special case m = n, i.e., where 
the scalar field energy density scales in the same way as 
the dominant component. The case where m = n is pro- 
duced by an exponential potential J0-p[, while negative 
power-law potentials give n < m H]^] ■ 

In this paper we provide a comprehensive classification 
of all solutions of this type, when the energy density is 
dominated by the perfect fluid. We show that the only 
potentials which lead to this sort of behavior are the 
previously studied exponential and negative power-law 
potentials, and a new class of positive power-law poten- 
tials. We study the general properties of such solutions, 
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including stability, and examine how well they might do 
in giving the desired cosmological behavior. 



II. SCALING SOLUTIONS 

A spatially-flat homogeneous Universe containing a 
perfect fluid with energy density p and pressure p, plus 
a scalar field (j> with potential energy V (</)), satisfies the 
equations 



H 2 



8ttG 



p=-3H(p + p), 



(1) 

(2) 



where dots are time derivatives. We will assume that the 
perfect fluid has equation of state p = (7 — l)p, which 
immediately implies 



p oc 



1 



37. 



The scalar field </> evolves according to 



-ZH<j)- 



dV 



The total scalar field energy is 

- v{<t>) + \i> 2 



(3) 



(4) 



(5) 



and we are interested in solutions for which p^ oc R n 
when p&_^L p and p oc R~ m . 

Eq. (||) allows the scalar field equation to be written 
in the useful form 



~3H<j) 2 



(6) 



If we divide Eq. (Q) by p^ and use p^jp^ = —n(R/R), 
then we obtain 



072 
P<t> 



n 
6 



(7) 



Thus, power-law behavior for the scalar field energy den- 
sity requires that the scalar field kinetic energy remain 
a fixed fraction of the total scalar field energy. The con- 
verse is true as well. This makes sense: if the kinetic 
energy evolves to become either dominant or negligible, 
then p^ will scale as 1/a 6 or remain constant, respec- 
tively. The former is not what we want, and the latter 
no different from a genuine cosmological constant. These 
two extreme cases also delimit the possible scaling behav- 
ior for the scalar field energy density: < n < 6, with the 
lower limit corresponding to potential energy domination 
and the upper limit to kinetic-energy domination. 



A. Exact Solutions 

Our basic method of solution is to assume the desired 
behavior of p^ and p and substitute into Eq. (|]) . A simi- 
lar procedure was first undertaken by Ratra and Peebles 
Q, who confined their attention to the cases of matter 
and radiation domination, m — 3, 4, and were interested 
in certain classes of solutions. Our development parallels 
and extends their analysis. 

When the perfect fluid with p oc R~ m dominates, 



R oc t 



2/m 



and Eq. (Eh becomes 
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dV 



(8) 



(9) 



The desired scaling behavior for p^, substituted into 
Eq. (0), gives 



4> oc t 



-n/m 



(10) 



Consider first the case m = n. Then Eq. (10) can 



be integrated to give cf> oc \n(t). Substituting this into 
Eq. (0) and solving for V((j>), we obtain 



(11) 



This is the well-investigated exponential potential ^ , for 
the limiting case where p^ ^ p. Although A can be 
positive or negative, those cases are physically identical, 
simply corresponding to a reflection of the </> trajectory 
about the vertical axis. 

Provided A 2 > m, the unique late-time attractor is a 
tracker solution with p^ = (m/\ 2 )p tota i |5|||. For exam- 
ple, the scalar field will rcdshift as 1/a 4 during radiation 
domination, and then switch to 1/a 3 once matter dom- 
ination commences. Although we derived it assuming 
P4> «C p, in fact this solution exists for any fractional 
scalar field density Q^, through the appropriate choice of 
A. 

While mathematically intriguing, such solutions seem 
uninteresting as candidates for a cosmological constant. 
During nucleosynthesis they behave as radiation and 
hence act like extra neutrino species, and are limited 
to < 0.2 during radiation domination and hence 
f2^ < 0.15 during matter domination, well below the 
desired density HtJ. A similar constraint arises from 
suppression of density perturbation growth (?j . Anyway, 
such a scalar field is presently evolving like matter and 
so won't explain the supernova measurements even if it 
were permitted with a more substantial density. 

These bounds can be evaded if the field does not enter 
the scaling regime until late in the cosmological evolution, 
e.g. after nucleosynthesis for the first bound, and after 
structure formation has been initiated for the second. 
However, this requires that the scalar field begin with 
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more or less its present density, and so provides no answer 
to the original fine-tuning problem. 

Now consider the case m ^ n. In this case, integrating 
Eq. (|l(3) yields 



<jy = At l - njm 



(12) 



The second integration constant has been absorbed by 
horizontal translation of <f>. Substituting the required be- 
havior into the scalar field equation leads to the potential 



V(</>) - A- 



( m) ( 
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2ri 



A 



where 



2/t 



(13) 



(14) 



The constant of integration, which would otherwise ap- 
pear in V{4>), vanishes because for scaling we need the 
kinetic energy to be a fraction n/6 of the scalar field en- 
ergy density. 

Scaling behavior can therefore occur provided the po- 
tential has a power-law form. If the exponent a is neg- 
ative, then m > n and the scalar field energy density 
grows compared to the matter, whereas if it is positive 
the opposite is true. We can rewrite Eq. ( O ) as 



a - 2 



in. 



(15) 



Since m and n are positive, Eq. (|15j) shows that scaling 
solutions exist for positive a only when a > 2 (in section 
III below, we consider what happens for a < 2). 

We have thus determined all potentials which give 
power-law scaling of when the dominant density com- 
ponent also scales as a power of R. The negative power- 
law and exponential potentials have been studied in de- 
tail B' our new resu lt i s the existence of scaling 
solutions with the positive power-law potentials. 

For most of these potentials, the differential equa- 
tion governing the evolution of (f>, Eq. ([)]), is nonlinear, 
and the solutions we have derived for <p(t) are particular 
rather than general solutions (in the study of nonlinear 
differential equations, these are known as 'singular solu- 
tions'). Hence, although there can be no other potentials 
which produce scaling behavior, there is as yet no guar- 
antee that the potentials we have derived produce general 
solutions (as opposed to singular solutions) which display 
the desired scaling behavior. Put another way, we must 
show that the singular solutions we have derived in this 
section are attractors of the equations of motion. 



B. Attractor structure 

The attractor structure of the exponential potential 
has been analyzed in detail elsewhere B so we will not 



concern ourselves with that potential here. The attractor 
structure of the negative power-law potentials has been 
discussed by Ratra and Peebles || for the cases m = 3,4. 
We extend their discussion to the case of arbitrary m, and 
also consider the case of positive power laws. 

We substitute a potential of the form V{4>) — Vo<p a into 
Eq. (||). However, note that the multiplicative constant 
in front of dV/d<j) can be absorbed into a rescaling of t. 
Henceforth, we assume such a rescaling and write 



6 l X 

-0 

m t 



(16) 



For this rescaled equation, the constant A in Eq. (TL2) is 



A = 



a-2 



6 
m 



a-2 



VO-2) 



(17) 



Note that A is well-defined, and the solution given by 
Eq. (O) valid, only for 
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ui a — 2 



6 — n 



>0, 



(18) 



which is satisfied automatically as long as n < 6. 

Following Ratra and Peebles 0, we make the change 
of variables 



t = e T ; u(t) 



Mr) ' 



(19) 



where 4>c(t) is the exact (singular) solution given by 
Eqs. and @. With these changes, Eq. @ be- 
comes 



4 6 
2 — a m 
2 



1 )u' + 



(20) 



a-2 



a — 2 m 



L ) = o, 



where the prime denotes the derivative with respect to 
t. This can be split into the autonomous system 



u =p, 



m 2 — a 



(21) 



a-2 



m a — 2 



For positive a, the interesting case is when a is an even 
integer, and then there are three critical points, all with 
p = and with u = —1, and 1. All three of these rep- 
resent solutions which asymptotically approach </> = 0, 
= 0. The u = +1 and u = —1 critical points, when 
they are attractors, represent solutions which asymptot- 
ically approach the singular solution. They give mirror- 
image trajectories; the +1 attractor represents solutions 
which go to 4> = from the positive <j) direction, while 
the —1 attractor gives solutions which approach = 
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from the negative <fi direction. The u = critical point 
corresponds to solutions in which <p goes to zero faster 
than in the exact solution in Eq. (p"2|). 

Linearizing Eqs. ( ^l|) about the u = 1, p = critical 
point and solving for the eigenvalues A± of small pertur- 
bations about this point, we find 



(22) 



a = 8 , m = 4 




2 - a 



a~2 



6 

77? 



For the cases m = 3 and m = 4, this equation reduces 
to the Ratra-Peebles results |j) . The behavior of these 
eigenvalues is somewhat clearer if written in terms of m 
and 7j, using Eq. (14): 



Ah 



2b — m — 6 ± ^/ (2n — 777 — 6) 2 + 8777,(7?, — 6) 



2m 



(23) 



The necessary and sufficient condition for stability is that 
the real part of both A+ and A_ be negative. If the quan- 
tity under the square root in Eq. (E3I) is negative, then 
this corresponds to the requirement that 2n — 777 — 6 < 
and gives a stable spiral. Note, however, that because 
?? < 6, the second term under the square root is al- 
ways negative. Hence, if the quantity under the square 
root is positive (so that both eigenvalues are real) then 
2n — m — 6 + y(2n — m — 6) 2 + 8777(77 — 6) < whenever 
2ti — 777 — 6 < 0. Hence, the condition for stability is 
just 2n — 777 — 6 < 0, regardless of the value of the quan- 
tity under the square root (although that will determine 
whether the stable singular point is a stable spiral or a 
stable node). In terms of a, the stability condition is 



a < 2 
a > 2 



6 + 777 
6 — 777 
6 + 777 
6 — 777 



Negative a ; 
Positive a . 



(24) 
(25) 



The first of these is always satisfied, showing that the 
scaling solution for the Ratra-Peebles potentials (a < 0) 
is a stable attractor for all values of a (as noted by Ratra 
and Peebles for 777 = 3,4). For positive a, however, the 
scaling solution is a stable attractor only for sufficiently 
large a. For example, in the matter-dominated era, at- 
tractor scaling solutions exist only for a > 6, while in the 
radiation-dominated era, this condition becomes a > 10. 
For 777 = 3, a = 6, and m = 4, a = 10, we have a vortex 
point at the singularity, which is neutrally stable but not 
an attractor. 




a = 10 , m = 4 




a = 12 , m = 4 




FIG. 1. Three phase planes for the radiation-dominated 
m = 4 case. From top to bottom, a = 8, 10, 12. The top 
panel spirals outwards, and the bottom one inwards. 



(a radiation-dominated universe). From top to bottom, 
C. Phase plane analysis the choices are for the singular solution [u = ±1, p = 0] 

to be unstable (a = 8), to be marginally stable (a = 10), 
A complete analysis requires a numerical solution. and to be stable (a = 12), respectively. To visualize the 
Three phase planes are shown in Fig. |for the case m = 4 physical meaning of these trajectories, note that trajec- 
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tories which cross the u = vertical axis correspond to 
solutions in which the <f> field is oscillating about the min- 
imum in the potential at (j> = 0, while trajectories which 
are confined to the left or right side of the u = axis 
correspond to solutions in which (j) goes to zero without 
oscillating. 

In the top figure, any point in phase space spirals out to 
infinity. This corresponds to a solution which never stops 
oscillating. The amplitude of the </> oscillations is decreas- 
ing; the spirals move outwards because they lose ampli- 
tude more slowly than the exact solution cf> c (r) . Although 
such a solution exhibits a form of scaling behavior (see 
section III.B below) the scaling exponent is not given by 
Eq. (15). The middle figure corresponds to marginal sta- 
bility and the phase trajectories are closed loops (which 
circulate clockwise), with the solution oscillating around 
the singular solution but not approaching it. The tra- 
jectories which cross the = axis represent solutions 
in which cj> oscillates forever. The trajectories closer to 
the singular point which do not cross the = axis rep- 
resent solutions in which does not oscillate about the 
minimum, but approaches zero as t x ~ n l m times an oscil- 
latory function. The bottom figure shows the attractor 
situation; depending on the initial conditions the trajec- 
tory may circulate several times around the two critical 
points (which corresponds to the field oscillating about 
the minimum) before circulating only about a single point 
(in which case the field stops oscillating and falls steadily 
toward the minimum). In the latter case the exact so- 
lution given by Eq. ( |l2| ) is multiplied by an oscillatory 
function which has a steadily decreasing amplitude. 

We can repeat our stability analysis for the u = 0, 
p = critical point. We find that the eigenvalues are 
both real, with A_ = (n — 6)/m, which is negative for 
n < 6, and A+ = n/m — 1, which is positive (since n > m 
for positive a). Thus, the (0, 0) critical point is always an 
unstable saddle point. This may seem bizarre, since the 
(0, 0) critical point corresponds to the field lying motion- 
less at the bottom of the potential. However, remember 
that the singular solutions also asymptotically reach the 
minimum, and our result simply means that if the field is 
perturbed slightly from this minimum, it returns to the 
minimum via the singular solution. 



pens, the Universe enters an inflationary regime, which 
has in fact been investigated in the early Universe con- 
text under the name "intermediate inflation" JTTj[ ] . The 
expansion rate asymptotically becomes 



i?ocexp i 4 /( 4 -«) 



(26) 



and the fluid becomes less and less relevant. The in- 
flationary regime may be preceded by a period of non- 
inflationary scalar field domination, if the scalar field 
comes to dominate while is sufficiently small. 

The scalar field density grows with respect to the fluid 
regardless of whether the Universe is radiation or mat- 
ter dominated, so these solutions do not exhibit a 'trig- 
gered' transition into the inflationary regime. Rather, 
the timing of that transition is governed by the initial 
conditions, and for the domination to be a recent event, 
one has to arrange for the initial scalar field density to 
be well below the radiation density. The tuning is not 
however as severe as with a pure cosmological constant, 
since the redshifting of the scalar field may be quite sim- 
ilar to that of the fluid ||. A particularly interesting 
case arises for a = —6; such a scalar field will scale 
as matter during the radiation-dominated era, and then 
grow relative to matter, as p^ oc i? -9 / 4 , once the matter- 
dominated era begins. If the field is generated initially 
with w Pmattcr *C Pradj then it will continue to evolve 
with w Pmattor until matter domination. The onset of 
matter domination then triggers a change in the evolu- 
tion of the scalar field energy density, and p^ begins to 
evolve in a manner close to a curvature density until it 
comes to dominate. 

An interesting question is whether it might be possible 
to find inflationary scenarios capable of providing suit- 
able initial conditions. A possible objection to the above 
is that in standard cosmological scenarios the energy den- 
sity which today is in non-relativistic particles (especially 
the baryons) starts out as highly relativistic, only later 
to change its equation of state on cooling, rather than 
already existing as a trace amount in the early Universe. 
However, the scenario just outlined bears some similarity 
to suggestions for creation of the cold dark matter at the 
end of the inflationary epoch . 



III. APPLICATIONS AND SPECIAL CASES 

A. Negative power-laws 

If a is negative we have a decaying power-law potential, 
in which the field can roll forever. These are the poten- 
tials first investigated by Ratra and Peebles (9j. They 
were recently reexamined in some detail by Zlatev et 
al. ||, in the context of the current observational situa- 
tion. Because the scalar field density grows relative to the 
fluid, eventually the approximation that the fluid energy 
density is dominant will break down. When that hap- 



B. Positive power- laws 

Although positive power-law potentials are more com- 
monly associated with driving an inflationary expansion, 
provided they satisfy Eq. ( |25|) then we have shown that 
they too permit stable scaling solutions. In this case 
to > n, and so the scalar field becomes progressively less 
important as the evolution proceeds, better and better 
justifying the neglect of the scalar field terms in the Fried- 
mann equation. The scaling solution for goes smoothly 
to zero as t — > oo, without oscillations. 

Note, however, that these potentials (for even a) can 
also support oscillatory behavior, with ]12j 
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Iog 10 R 

FIG. 2. The evolution of the energy densities in a radia- 
tion-dominated Universe with a — 30. The vertical axis is in 
arbitrary units. The dotted line is the radiation energy den- 
sity, and the solid line is the scalar field potential energy. The 
scalar field kinetic energy, shown as the dashed line, oscillates 
through zero out of phase with the potential energy, but does 
so too sharply for the plotting resolution. 



oc R 



-6a/(a+2) 



(27) 



Despite the power-law behavior, these solutions are not 
encompassed in our definition of scaling, as the scaling 
law arises only after averaging over oscillations, while 
within each oscillation, energy is continually being con- 
verted between potential and kinetic. When the field is 
oscillating, and Eq. ( p7|) applies, the scaling of p<p with R 
is independent of the equation of state of the dominant 
component of the density; this differs from the attractor 
solution in which the scaling of p^, with R depends on 
77i. Furthermore, such solutions apply to oscillating fields 
even when the scalar field density itself is dominant. In 
our phase diagram, Fig. [j], this oscillating solution corre- 
sponds to the regime in which the phase space trajectory 
winds around both attractors. 

Whether the oscillatory behavior perseveres depends 
on the stability condition of Eq. (25), which if satisfied 
implies that the scaling behavior found in Section 2 is 
the attractor. The oscillating solution has an amplitude 
0max ~ R~ 6 /( a+2 \ which matches the redshift depen- 
dence of the singular scaling solution if the stability con- 
dition is saturated, leading to the closed loops seen in the 
middle panel of Fig. |l|. 

If the stability condition is satisfied, this can lead to 
some interesting behavior. Consider the limit where a is 
very large and positive in the radiation-dominated era. 
If initially p^, p ra d, but 4> is oscillating rapidly, then 
the <j) energy density will scale roughly as p$ oc l/R 6 



and eventually fall below the radiation density. When 
this happens, however, the scaling behavior will take over 
and the cf> energy density will scale as 1/ i? 4+e with K<1. 
We show this evolution with a = 30 in Fig. ||. Initially 
we have inflation, and then the field undergoes oscilla- 
tions which are heavily dominated by the kinetic energy, 
so its energy falls off at nearly l/R 6 , similar to kina- 
tion [jlj|. Finally the scalar field becomes subdominant 
and stops oscillating, instead entering the scaling solu- 
tion with p^ oc 1/ 'i?( 4+2 / 7 ). This scenario provides yet 
another "natural" mechanism to give a scaling solution 
with roughly equal to the density of the dominant 
component, since the scalar field density drops rapidly 
relative to the radiation density until p^ « p ra( j, after 
which p0 decreases at nearly the same rate as p ra d- A 
similar behavior has also been noted for the case of neg- 
ative power-laws when the <p density is initially much 
larger than its attractor value the difference in the 
case discussed here is that the attractor is reached even 
when p0 /7 ra d initially. 

If the stability condition Eq. (^J) is n °t satisfied, then 
the oscillations continue indefinitely as in the top panel of 
Fig. |l|. Depending on parameters, the scalar field energy 
density may be either increasing or decreasing relative 
to the fluid energy density, e.g. for the choices in the 
figure, p ra( j oc l/R 4 while p^, oc l/R , so the scalar field 
becomes less and less important. 

What happens for < a < 21 For a = 1 or 2 it is easy 
to find the exact solutions. First consider a = 1. Then 
Eq. (M) becomes 



7 6 1; 
777, t 

and the exact solution is 

6 = A + Bt 1 - 6 /" 1 



1 = 0. 



1 m 



2 m + 6 



(28) 



(29) 



where A and B are constants to be determined by initial 
conditions. This solution does not display scaling behav- 
ior, and, not surprisingly, 4> — > — oo as t — > oc, so it is of 
little physical interest. 

The m = 2 case is more interesting. For this case, we 
get the linear equation 



6 1 

777 t 



cp + <f> = 0. 



Taking cf> = Ot 1 ^ 2 3 / m , this equation reduces to 



t 2 9 + te + t z [1 - (1/2 -3/m) 2 ] 6 = 0. 



(30) 



(31) 



This is Bessel's equation of order 1/2 — 3/m\, so the 
general solution for <f> is 



= t (l/2-3/m) [ A J u {t) + BN v {t)] 



(32) 



where A and B are constants determined by the initial 
conditions, and J„ and N v are Bessel functions of order 
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is, with v = 1 1/2 — 3/m|. In the limit of large i, the solu- 
tions in Eq. (|32|) all oscillate sinusoidally, with amplitude 
decaying as i~ 3 / m , so p^ oc R~ 3 . Thus, the solutions for 
this potential always oscillate, and the density scales as 
in Eq. ©. 

Finally, some exact solutions exist for the case to = 3 
(matter domination), for which Eq. ( [l6| ) reduces to the 
Lane-Emden equation. This equation can be solved ex- 
actly |T^] for a = 1,2 and 6, with the Lane-Emden 
boundary conditions corresponding to the initial con- 
dition 4> = at t = 0. The solutions for a = 1 and 
a = 2 are special cases of the solutions discussed above. 
The case a = 6 corresponds to the transition between 
attractor and non-attractor behavior and represents the 
analogue (for the matter-dominated case) of the a = 10 
potential for the radiation-dominated universe shown in 
Fig. 0. The phase diagram in this case resembles the 
middle diagram in Fig. |l|. However, the Lane-Emden 
solution is not the u = ±1 attractor solution; rather it 
corresponds to the unstable singular point at u = 0. This 
arises because the Lane-Emden boundary conditions cor- 
respond to initial conditions which lie exactly on the 
singular point u — 0, p = 0, and the solution remains 
there as t — > oo. The Lane-Emden boundary conditions 
are unphysical when applied to the scalar field evolution 
equation, since t — is undefined in the cosmological 
context. 



C. The ZWS potential 

Zlatev et al. || made a detailed analysis of the rather 
unusual potential 

V{4>) oc exp fe) - 1 . (33) 

This potential is introduced in recognition of the fact 
that the simple power laws do not exhibit the ideal cos- 
mological behavior, in that the scalar field density grows 
relative to matter during the matter-dominated era only 
if it also grows relative to radiation during the radiation- 
dominated era. While the fine-tuning problem of why 
the cosmological constant took so long to dominate is 
certainly less severe with these power-law potentials than 
with a pure A term, it still remains and one requires ei- 
ther an extremely low density in the scalar field at early 
times, or alternatively to have n extremely close to to 
(i.e. very large positive or negative a) so that the scalar 
field requires a very long time to catch up with the con- 
ventional matter, say from an initial state of equiparti- 
tion with a large number of fluid components. Further, 
the latter resolution, while superficially attractive, will 
fail for the same reasons that the exponential potential 
does, namely nucleosynthesis and structure formation; its 
dynamics are extremely close to the exponential case. 

The purpose of a potential such as in Eq. ( p3| ) is to 
change the slope of the scalar field, and hence alter 



the character of the scaling solutions with epoch. For 
mpi/</> <C 1, the potential decreases more rapidly than 
any power law. This initial steepness guarantees n ~ to, 
and the field is drawn to this approximate tracker behav- 
ior. Later, when (j) 3> mpi, this potential asymptotically 
approaches the form V{<jj) oc 1/0, with the scaling solu- 
tion n = m/3, producing scalar field domination. This 
strategy certainly does yield the attractive observational 
consequences explored by Zlatev et al. ||. The draw- 
back is that the change of behavior is now governed by 
the form of the potential, and not primarily by the equa- 
tion of state of the accompanying fluid. That the scalar 
field begins to change its behavior around the epoch of 
matter-radiation equality is because the feature of chang- 
ing steepness in the potential has been placed in the ap- 
propriate place. This represents tuning of a different sort 
to the usual tuning of A models, but a tuning nonethe- 
less. 



IV. DISCUSSION 

Our results indicate that exact solutions for the scalar 
field, which give scaling behavior when the expansion of 
the universe is driven by a dominant component with 
density pdominant, are possible for only three classes of 
potentials: 

1. Exponential potentials (p^ scales as ^dominant)- 

2. Negative power-law potentials (p<^ decreases less 
rapidly than ^dominant ) ■ 

3. Positive power-law potentials (p^ decreases more 
rapidly than ^dominant)- 

The first two cases have been extensively discussed else- 
where; the existence of the third class is our major 
new result. The negative power-law potentials V oc <jf 
have attractor solutions for all values of the exponent a, 
while the positive power-law solutions require an expo- 
nent a > 2(6 -I- to)/ (6 — to) for attractor behavior to 
occur. 

Our results do have one practical limitation: we have 
confined our attention to exact solutions. It is cer- 
tainly possible, for example, for approximate solutions 
to exist which are very close to scaling behavior, e.g., 
Pc/, oc R~ n f(R), where f{R) is a slowly-varying func- 
tion of R. If f(R) varies sufficiently slowly, then there 
may be no practical distinction between a solution of this 
type and our exact solutions. The particular potential of 
Zlatev et al. is of this type. It is not practical to system- 
atically classify all such approximate scaling solutions, al- 
though it may be possible to provide conditions on V(4>) 
which allow for such solutions Eal. 
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